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SUMMARY 

bp«.n condition potential flow calculation method has 

been combined with a momentum Integr-il boundary layer method and a basrflow 

airSey^nd's'JSf vlscous-lnilscld Inte^er^ncr^p 
seoaraJlM «nd Jc ’resultant program considers laminar and turbulent 
eparatlon and Is, therefore, applicable to thin or thick airfoil stall 

the sepa^tJoi lu^le 

recombination region behind the airfoil. 

th« y*"® distributions and equivalent airfoil shapes. Including 

thickness of the viscous regions, are compared with flow 

fnd'waK"L“nt:SJLf°^ dlsp^acemenrth "nesses 

suLests the use calculated shape. The comparison also 

suggests the use of the analytical solution to evaluate the measurements. 


INTRODUCTION 


of alKrafrafhl.r!S?L'°J °f t*-* aerodynamic characteristics 

aircraft at high angles of attack, the problem of determining the flow field 

around a stalled wing was addressed. A necessary prellmlnarrsLp Ls bLn tJe 

M veru, I^^aSn^c^MfL f developed primarily 

extenainn f bhe theoretical components prior to their 

extension to a three dimensional calculation method. 

Three types of boundary layer separation have been Identified as centri- 
cal ”70^0! stall in combination. The classl- 

laver can he e<rh ^1 trailing edge separation. The separating boundary 

a JSil JralllifeJ d®' turbulent. Separation first occurs at the 
airfoil trailing edge and moves forward with Increasing angle of attack ‘ This 

(t/c > isf^mr f relatively thick airfoils ' 

fiLl L 7 form of separation Is short bubble laminar separa- 

tirbuTenrf} J separation Is followed almost Immediately by transition to 
turbulent flow and boundary layer reattachment. This bubble develops as ^h^ 

uiS JeatJach^fn^n®*'? increased, or as Reynolds number is decreased, 

^ longer occurs. The bubble then is said to have "burst” 
n thick airfoils short bubble laminar separation causes a thicker downstream* 

diJ^cJly JHlrfoll Stan!"”®’' airfoils, short bubble bursting can lead 
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with 't^Tnol Jowe^suSL! b^M fT recombines 

under certain circumstances of airfoll^thlrV^ * trailing edge. However, 

reattaches to the airfoirLrLce fn^ " ^ Reynolds number, the flow 

This is referred trls s lnnl ff ^ ”®'' turbulent boundary layer is formed, 

attack cause the bubble lentfth^^*?«’^ separation bubble. Increases of angle of 
airfoil trailing edge. * increase until reattachment moves behind the 

upon fpo?e"ntiaJ progression is based 

or pressure distribution ^o Ka a *" either the airfoil shape 

Of the a^foil s^face ms ooLl.f “3r portL 
with momentum Integral ‘bo^darJ i ^ ^ calculation method is combined 
analysis da«rlKlS of “l-=“laMon "athods and a conponant 

Interactlortteow ?ha ‘d for. a onifled yiscous-lnvlscld 

included uhether they occur sln8iror°ln°cOTblM«oor “"‘’''d “e 


NOMENCLATURE 

airfoil chord 
pressure coefficient 
airfoil thickness 
velocity on airfoil surface 
freestream velocity 
horizontal distance 
angle of attack 


MATHEMATICAL MODEL 


..odelId%rfn;;8”rtrecul«UnrdJLr°^^ " “paratlon is 

around the airtoU ™d l«o ““ 

^n^ r -- - --- tE:ir 

placement thickness Is added to the beslf eUfoll ?o f ■"<' d‘a- 

pressure predicted by aTas. ?^^^^^^^ thickness is not known, so the 

The output of this calculation Is m “ '*'! '’"“"‘‘•ty condition, 

lent airfoil shape. The updated' pressSre dls?r“n"hen1rired“f tSe“‘''°' 


Input for Inproved viscous flow calculations. This proce.iute Is raoeatad 

tiv~E£S 


Potential Flow Calculation 

recognized that when the boundary conditions of a 
calculation method a7e linearized ?n the 


Separation Bubble Model 

.a. T-r. r (^.rEa'nSir" 

a“„d“t^^^Zarra““r^:l:tfnf ° 

OII^OT- fi s-u^ limiting streamline between the Inner and 

o? tL '•■a °utar flow Is assured? iSa 

babble ilTC? i? determines the pressure rise across the separation 

or recombination °' "attachment 

to increase towaJd'?h%“‘?|s?ut a^“^e\t^'X"n^^1^e^“X'>ati'^n^ T?s“«?, 

iT-r ?x?r ?hf - - bubM? 

is sh?m“?“n8?r?'2?"'l^?hls“?L?'’?J''’' w“ "**' aaParation 

s'?T?‘S!? T‘"l T ''>a^?w«'su????e'mS f??m.“'’?"ak?'''“he 

point new turbulent boundary layer is formed which continues dowilstream. 
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F'‘om this physical picture, it can be seen that the pressure within the 
separation bubble is limited by the static pressure at the closure of the 
bubble and not by the pressure at separation. A solution to the flow problem 
is reached when the presence of the separation bubble modifies the pressure 
distribution around the airfoil sufficiently that the pressure at separation 
agrees with the bubble pressure. 


VERIFICATION OF THE METHOD 

The Generalized Boundary Condition potential flow calculation method and 
the separation bubble model were combined with momentum integral boundary 
layer methods into a procedure for calculating the viscous-lnviscld inter- 
actions for a wide range of airfoils. The resulting program is self-contained 
and includes all of the logic for distinguishing between the types or combina- 
tions of separated flow that are present. An iterative procedure is used to 
go from the Initial inviscid solution to the viscous solution. In addition, 
iterative sub-loops are required to determine the proper bubble pressures for 
the long laminar separation bubble and the trailing edge bubble. The iteration 
procedure is described in Reference 1 . 

In order to test the ability of this analytical model to predict the flow 
field around an airfoil with various types of separated flow present, five 
airfoils of different thickness ratios were studied. The experimental data 
for these airfoils were taken from References 5 and 6 . They^were; ( 1 ) the 
633-018 airfoil with stall development by progression of turbulent separation 
from the trailing edge forward, ( 2 ) the 63i-012 airfoil which stalled by the 
sudden bursting of a short separation bubble, ( 3 ) the 63-009 airfoil that also 
stalled by short bubble bursting, (4) the 64A-006 airfoil with stall develop- 
ment by the progressive growth of a long separation bubble until its reattach- 
ment point moved past the airfoil trailing edge, and (5) the GA(W )-1 airfoil 
which stalled similarly to the 633 -OI 8 airfoil but represented more advanced 
airfoil design. The calculations were made for the specific angles of attack 
at which the surface pressures on the model had been measured. 

Representative calculated pressure distributions and equivalent airfoil 
surfaces are shown in Figures 3 and 4. The inviscid pressure distribution is 
shown for comparison with the calculated pressure distribution and the measured 
pressures. The equivalent airfoil shape is the combination of the airfoil and 
the displacement surface of the viscous flow. The limiting streamline between 
the continuing viscous flow and the recirculating bubble flow also is shown. 

Figure 3 shows a thin airfoil with a long laminar separation bubble. In 
order to simulate this flow, it was necessary to assume a bubble pressure and 
a point of reattachment. The pressure distribution was prescribed between 
separation and reattachment, and the airfoil shape was prescribed outside of 
this region. The subsequent equivalent airfoil shape was determined by the 
Generalized Boundary Condition program and the shape of the limiting stream- 
line from the base flow theory. Reattachment was defined as the point where 
the Umlting streamline intersected the airfoil surface. If the pressure rise 
between separation and reattachment corresponded to the predictions of the 
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pressure was assumed, ai^the^asi^erreaTtTc^^^^^^ ^ d 

results of the previous calculation. was updated by the 

The predicted bubble pressures shown -in Ploure av^ u 
experiment and there is a k<nk <n ^^ 8 ^® 3 are somewhat lower than 

Ih. first effect ‘rdsa L tS^ ?*=«“’"“»" « resttechn»„t. 

Shear layer ,rcth across the separstlon biMe'u :SjJaXi"?ert! 

Figure 4 shows the calculated results for the NACA 6 ln ni 9 -f ^ 4 i 
angle of attack that is sufficientlv u i o3x-012 airfoil at an 

bursting has occurred, Thirfiau^e laminar separation bubble 

the equivalent alrSli c^tLr °f the fact that 

therefore has the pot^M^l f !? ?" “thod 

the lift and pltchf™ iSis uaf "f “ 

not considered to he*sufflclmtw*!c tlote the present method Is 

However, the mettod a sSrfLr t f ‘7 *“'* “ parameter, 

capable of this calculation. ' hb“ 8 slarlty potential flow method should be 

The equivalent alrfoU surfaces of the NASA rA^-m i 41 

lated at three angles of attack and are shown “ Flgwes 5 6 LdT 

Imposed on the flow field measurements of Seetharam aJd Wents’frrJ 
measured boundary laver dlqnl 9 f><aman»- e-im^ i (ref. 6 ). The 

shown. It was tL obiecrof thir^Jf centerlines also are 

just how the bubble matchpH th develop some insight on 

for the development of a surface ^slnBularltv^m'^rh d^^r Importance 

advantage of allowing an evaluation of rh*> ^ should have the further 

it Should give an app%e“i:«rot:h:'e?J:crof”S^ ^ 

tlon h"bh”e" haf^« Jet fJjLf T °5 thPaf- 

».y to the tr"u“ Jdje H»erer th«e' 1 s L“* 

the trailing «dge. The measureT7lsVl“ e^ent th^^^^^ 

the equivalent airfoil surface Hnpc u agrees reasonably well with 

wake is much thicker thjrthj m^.. ; a ! "f ! ibe calculated 

did not allow sufjLuSrpJJjsJJeTl^eSr^herkr*"' 

(Flgu« t) separation bubble has formed 

with the calculated eJJlJalJnt JlJtJll'* thickness does not agree 

centerline does. It if noted rh^ surf ace. although the meusured wake 

do not agree with the measured separaVlOT po“l 5 it'^ ‘*mis“'iT s”t' thicknesses 

tL^1\\fd&JnTr& 
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A fully developed trailing edge separation bubble at an angle of attack 
a » 1844 is shown as Figure 7# Of particular note is the constant pressure 
region that brackets the equivalent airfoil contour# The measured displacement 
thickness also agrees with the calculated contour. However, there is a marked 
deviation of the measured wake centerline. This Is probably due to the effects 
of the wind tunnel floor and celling which are known to have a pronounced 
on the airfoil downwash when separated flow is present. 

These results suggest an approach to the determination of wind tunnel wall 
corrections for separated ilows. Potential flow calculation methods have been 
developed to simulate the flow around airfoils between parallel walls. The 
addition of singularity panels simulating such parallel walls in the current pro- 
gram should lead to closer agreement with the measured flow shown in Figure 7. 

Th® difference between the two cases would yield both streamline curvature 
and blockage corrections. 


CONCLUDING REMARKS 


The McDonnell Aircraft Company Generalized Boundary Condition potential 
flow calculation method has been combined with momentum integral boundary 
layer methods and a component analysis base flow theory to develop a method 
for predicting vlscous-lnviscld Interacting flows on airfoils beyond the 
appearance of boundary layer separation. The physical models of the several 
phases of such flows have been Identified and combined into a functioning 
whole that accounts for the Interaction of attached flow, short and long 
leading edge separation bubbles, and trailing edge bubbles. The resultant 
computer program has been used to calculate the pressure distribution and 
equivalent airfoil shape for airfoils exhibiting the different types of 
separated flow. It also has been used to evaluate flow field measurements 
around an airfoil with trailing edge separation. 
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Figure 1.- Applications of generalized boundary conditions to separated 

flow modeling. 



Figure 2.- Model of boundary-layer separation bubble developed from the 
component -ana lysis base-flow theory. 
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CQUIVAUNT COMTOtiR 


Figure 3 


Figure 4. 



Prediction of airfoil pressure distribution and equivalent 
contour 64A-006 airfoil, a « 7.31^*. 



- Prediction of airfoil pressure distribution and equivalent 
contour 63,-012 airfoil, a ■ 
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Figure 6.- Superposition of calculated equivalent airfoil shape on measured 
flow field of GA(W)-1 airfoil at a ■ 14.4®. 





Figure 7.- Superposition of calculated equivalent airfoil shape on measured 
flow field of GA(W)-1 airfoil at a - 18.4*. 



